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Abstract 

We put forward an improved version of the Galilean Genesis model that addresses the problem of super- 
luminality. We demote the full conformal group to Poincare symmetry plus dilations, supplemented with 
approximate galilean shift invariance in the UV and at small field values. In this way fluctuations around 
the NEC- violating cosmological background are made substantially subluminal, and superluminality cannot 
be reached by any small change of the solution, in contrast with the original model. Dilation invariance still 
protects the scale-invariance of correlation functions of a massless test scalar — which is the source of the 
observed cosmological fluctuations — but the explicit breaking of the conformal group can be potentially 
observed in higher-order correlators. We also highlight a subtlety in matching the NEC-violating phase with 
the standard cosmological evolution, and discuss the possible couplings of the Galileon to gravity. 



1 Introduction 

There is a recurrent yet somewhat vague connection between the Null Energy Condition (NEC) for the stress- 
energy tensor and the subluminality of signal propagation, (i) The NEC is a close relative of the dominant 
energy condition, which states that for any observer the energy-momentum density should be a time-like 
future directed vector — that is, no energy should leave the observer's future lightcone. (ii) In general 
relativity the NEC implies that no closed time-like curve can form, provided suitable boundary conditions 
are satisfied [1]. (Hi) Perhaps related to this, for linearized gravity in asymptotically flat spacetime, the 
NEC ensures that the light-cone defined by the gravitational field is always narrower than that defined by 
the underlying Minkowski spacetime [2, 3]. (iv) For effective theories involving scalar fields only, at lowest 
derivative level the subluminality of excitations about a given solution implies that the solution obeys the 
NEC [4]. (v) For the conformal Galileon [5], certain solutions violate the NEC, others feature superluminal 
excitations [6, 7]. In each of these instances the connection is far from vague, yet it is of a different nature 
every time. Moreover, it is very indirect in the example of (v), where the NEC- violating solutions are 
perfectly sensible yet the same effective theory admits other solutions (including small perturbations of the 
NEC-violating ones) that unavoidably exhibit superluminality. These facts hint at either of two opposite 
conclusions: (1) That the NEC/subluminality connection is generic and much deeper than shown by each 
individual example (and so far we have not been able to make it explicit); (2) That it is accidental, and 
peculiar to the examples above. 
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In this paper, we sever such a connection. We wiU exhibit an effective field theory (EFT) with the 
following properties: 

1. There exists a homogenous and isotropic NEC-violating solution that is stable against small pertur- 
bations. 

2. About this solution, perturbations are generously sub-luminal, with (5c/c ^ 1. This makes their 
subluminality particularly robust — it will survive generic deformations of the background solution, 
as long as these are not too large. This is in sharp contrast with the NEC- violating conformal Galileon 
of [6, 7]. 

3. The structure of our effective Lagrangian is protected against large quantum corrections by approxi- 
mate symmetries. 

Ideally, one would like to supplement these properties with: 

4. A Poincare invariant solution, also stable against small perturbations. 

5. A Lorentz invariant S'-matrix about this solution, obeying the standard positivity constraints coming 
from relativistic dispersion relations. 

6. Robust subluminality about this solution, in contrast with the NEC- violating conformal Galileon [6, 7]. 

That is, one would like to have a perfectly well behaved relativistic field theory that can be defined starting 
from a Poincare invariant vacuum state, and that features NEC-violating solutions that are also perfectly well 
behaved. Unfortunately, at least in our simple framework, this will turn out not to be possible: conditions 
4-6 contradict conditions 1-3. Notice however that the absence of a Poincare invariant vacuum state is 
not necessarily an inconsistency for a relativistic effective field theory — and we have obvious empirical 
evidence for this: there are systems like solids and fluids, which (a) exist as Lorentz-breaking states in 
a Lorentz- invariant microscopic theory (the standard model of particle physics); (b) can be described as 
Lorentz-breaking classical solutions in certain relativistic effective field theories [8, 9], which however cannot 
be consistently extrapolated to zero density and used to describe the Poincare invariant vacuum of the 
microscopic theory. For instance, we cannot hope that the vacuum of the standard model be well described 
by hydrodynamical equations. 

Our starting point will be the conformal Galileon [6] . We will build upon the results and the analyses of 
[6, 7], to which we refer the reader for details about the original model and for the notation. There, the reason 
we could not avoid superluminality was essentially one of symmetry. The high-degree of symmetry of the 
action — 50(4, 2) — is partially broken by the NEC- violating solution, leaving the de Sitter group S'0(4, 1) 
as the residual symmetry group. This degree of residual symmetry is enough to guarantee that in the UV — 
at scales much shorter than the 'curvature' scale of this de Sitter solution — excitations travel exactly at the 
speed of light, i.e., on the verge of superluminality. Then, it is a matter of details whether small deformations 
of the background solution will admit superluminal excitations. For the conformal Galileon with nonzero C3, 
this is always the case, and one can check that the level of superluminality is detectable within the effective 
theory [6] . Setting C3 to zero does not seem to improve the situation (we analyze this possibility at the end of 
Appendix C). One is thus led to consider the possibility of making the excitations of the background solution 
strictly subluminal, which can only be achieved by trimming the residual symmetries. One could consider less 
symmetric background solutions, which in our case would complicate many computations, or less symmetric 
actions, which is the direction that we will take. Obvious symmetries to jettison are the special conformal 
transformations, since the remainder — 4D Poincare and dilations — close into a subgroup of 50(4, 2). In 
Sees. 2 and 3 we will show that the resulting theory obeys properties 1-3 above. 

Once minimally coupled to gravity, our system can drive the Galilean Genesis phase of [7] . As far as the 
background evolution and the power spectrum of perturbations are concerned, the cosmological implications 
of our new model are the same as for the original one. However, the models become observationally distin- 
guishable from each other and from inflation at the level of higher-point correlation functions. We discuss 
the symmetry reasons behind these facts in Sec. 4. 

In Sec. 5 we discuss an unexpected subtlety: in the matching of our Galilean Genesis phase to a standard 
radiation-dominated one — that is, at reheating — not all the energy built up by violating the NEC can be 
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passed on to the subsequent phases. Even for an instantaneous reheating, as a consequence of the higher- 
derivative structure of our interactions, an order-one fraction of the energy gets lost — without ever violating 
stress-energy conservation of course: the loss of energy is due to a singularity in the equation of state at 
reheating, which induces a sudden, finite redshift of the energy density. In practice this does not change 
the fact that we can start with vanishing energy in Minkowski space, and generate an expanding universe 
full of energy whose post-genesis cosmology is the same as our universe's. However it introduces a novel 
question for models that, like ours, use NEC-violations as alternatives to inflation: whether the energy and 
the expansion rate created while violating the NEC can be inherited by the standard cosmology that comes 
after. 

In Sec. 6 we turn to generalizations of our scenario, involving higher order Galileons and non-minimal 
couplings with gravity, like for instance those of the 'covariant' Galileon of [10]. We argue that all these 
choices are radiatively stable, and we discuss in what circumstances some of the non-minimal couplings 
might be improvements of the effective theory. An exhaustive analysis of all these generalizations is beyond 
the scope of our paper, but for the cases that we can analyze straightforwardly, we argue that nothing is 
gained by considering these more general possibilities. A number of technical derivations are collected in the 
Appendices. 

To conclude our introductory remarks, we should also mention the ghost condensate case [11], and explain 
in what sense we are improving on it. There, one can have consistent NEC violations without superluminality 
[12] — in fact excitations are always extremely sub-luminal. Notice that, like in our case, the theory does 
not admit a Poincare invariant vacuum within the same effective field theory as the ghost condensate point ^. 
However, in the absence of gravity, the NEC-violating branch of the ghost condensate is unstable on large 
scales. This might be irrelevant in practice for cosmological applications, as long as the instability scale is 
larger than the Hubble radius, but it is certainly a substantial difference with our case, which features no 
instability whatsoever, even in the absence of gravity. 

2 The NEC-violating solution 

We consider a deformation of the original Genesis Lagrangian [7] 



(1) 



with A <C / and a a new dimensionless parameter of order unity^. We neglect gravity for the moment. For 
a = we recover the minimal Lagrangian studied in [7] which non-linear ly realizes full conformal invariance. 
For a conformal invariance is explicitly broken, but dilation invariance is preserved. We will address 
the radiative stability of this structure in Section 3. For the moment, we notice that for small tt the theory 
reduces to the ordinary cubic Galileon of [13, 14], and so dilation invariance is enhanced to internal Galilean 
invariance, tt — >■ tt -f b^^x^ . This statement is radiatively stable as the renormalization of the operators is 
local in field space, so that the action will remain Galilean invariant near the limit, and it will enjoy there 
all the standard non-renormalization properties of Galilean operators [13, 15, 16]: in particular, quantum 
corrections will not generate the higher order galilean operators, nor will they renormalize the coefficients of 
{dirf and (dTrfUTT. 

We are after a 'de Sitter' solution, where e^'^ takes the form of the conformal factor for de Sitter 
space [6, 7, 17]: 



Hot ' 

with a constant. For the Lagrangian (1) we find 



-oo < t < , (2) 



^2 _ 2 1 A^ 



^Both classical solutions might be allowed, formally, but they are always separated in field space by a region with 
jhosts, signaling the breakdown of the effective theory. 

^We use the mostly plus signature, so the quadratic term in (1) has the ghostly sign. 
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The gravitational stress-energy tensor associated with the action (1), computed from the action as 
^ '^'^77 with minimal coupling, is that of [7], 



T — — 



g Sgi' 

conformal 



t: 



11 
A3 

2A3 



[4{d7T)^d^TTd,7T~g^,{d7Tf] , 



supplemented with 



(4) 



(5) 



As we discuss in Section 6 and Appendix A, there are ambiguities in how the higher-order Galileons are to 
be coupled to gravity which can affect the stress tensor, but these ambiguities do not enter here. 
For a time-dependent solution the total energy density and pressure are 



P 
P 



3(1 + «) / 4 

e TT —t: 



A3 



e rrrVr H — — TT 

2 A3 3 A3 d< ^ 



(6) 
(7) 



where we neglected the effect of gravity on the background solution, a good approximation at early times. 

A time-dependent solution violates the NEC if and only if the combination p + p is negative. For the de 
Sitter solution (3) we find 

2/2 3 + a 



which is negative for3 



a > — 1 or a < —3 



(NEC violation) 



(8) 



(9) 



Expanding the action to quadratic order about the de Sitter solution (3), we get the free action for the 
perturbations ip = tt — vTds , 



quad 



3(1 + a) 



(10) 



The overall coefficient is positive-definite. However, to ban instabilities, we have to make sure that the 
coefficient in front of the gradient energy, the signal propagation speed squared. 



3(1 + a) 



is also positive. This is the case for 



1< a < 3 



(stability). 



(11) 



(12) 



Finally — and this is the novelty with respect to [7] — we can consistently require that perturbations 
propagate strictly subluminally, < 1. This happens for 



a > (subluminality). 
Demanding (9), (12) and (13), we find that for 

< a < 3 (NEC violating, stable and subluminal). 



(13) 



(14) 



^Note that for the conformal case a = 0, the sign of p + p depends only on the sign of the kinetic term, and is not 
affected by the cubic Galileon terms. As we discuss in Appendix A, this property persists for the higher Galileons 
with the right gravitational couplings. 
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the system violates the NEC, is stable against small perturbations, and these propagate at subluminal speeds. 
For much larger values of a (at fixed / and -ffo); the cubic term becomes irrelevant, as it is the only one 
which does not scale with a on the solution; in this limit the violation of the NEC must be associated 
with an instability [4], and, indeed, the system exhibits a gradient instability. If a is 0(1) the speed of tt 
excitations is substantially less than unity, and, by continuity, generic perturbations of the de Sitter solution 
— including non-infinitesimal ones — cannot spoil this property and make the lightcone superluminal. 

Notice that in our quick stability analysis above we have neglected the mass term in (10). Its precise 
value is enforced by non-linearly realized time-translational invariance [6, 17, 18], and it appears to have the 
wrong (i.e., unstable) sign. Nevertheless, as shown in [6, 19] the mass term in (10) is unambiguously not 
associated with an instability. The growing mode solution represents a constant time shift of the background, 
indicating that the background solution is in fact a dynamical attractor. More generally, whenever the size 
of the mass term is the same as the time-variation rate of the background solution, the 'mass' of excitations 
is ill-defined. For instance, we can make the mass term vanish via the field redefinition ^ = if/TTds — shifts 
by a constant under the spontaneously-broken time-translational invariance, which then forbids any mass 
term for it — or we can flip its sign via some other time-dependent field-redefinition. This is related to the 
fact that an instability that develops on a time scale that is comparable to that over which the background 
solution itself changes by order one, cannot be unambiguously called an instability. 



3 Radiative stability 

We now discuss in more detail the radiative stability of the action (1) and of the NEC-violating solution. 
The action is invariant under dilations, so that, neglecting gravity, only dilation-invariant operators will be 
generated. The most generic dilation-invariant operator is schematically of the form 

0,„„=.e(4-™-")-=/^'?^^|^, (15) 

where m and n are (non-negative) integers, tTc is the canonically normalized field, tTc = /tt, and the powers 
of A are such as to give the operator overall mass-dimension four. The dimensionless couplings with which 
operators of this form appear in the classical Lagrangian or get generated at quantum level are not necessarily 
of order one — like the four derivative quartic term in (1), they can be suppressed or enhanced by suitable 
powers of A//. Let us take this possibility explicitly into account, by defining the operators 

where q can in principle have either sign. 

We now want to understand whether the operators that get induced quantum mechanically modify 
substantially the dynamics inferred from the classical Lagrangian (1). A way to assess this, is to estimate 
their size on the classical solution (2), (3), and to compare it to that of the classical Lagrangian terms of (1). 
If we evaluate a generic Om,n,q on the solution we get 

The scaling is a consequence of dilation invariance. In the above notation, the classical Lagrangian (1) 
reads 

/:-Oo,2,0 + Ol,3.0+Oo,4,l , (18) 

SO that all terms have k = 0, and are all comparable'^ on the solution, scaling as 1/i^ • (//A)^. We now 
want to show that the only operators that get generated by loop diagrams have A: > 1, and are therefore a 
negligible correction to the classical Lagrangian. 



*We are here interested in the typical size of each operator and we thus disregard that the contributions to the n 
equation of motion of the operator (c'7r)^n7r cancel on our NEC- violating solution above. 
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Figure 1: Generic 1-PI loop diagram. The interaction vertices can he extracted from eq. (1), upon expanding 
in IT. 



Consider a generic IPI graph, like for instance that depicted in Fig. 1. The interaction vertices are those 
of (1), after we expand the exponential in powers of tTc- For any given graph, vertices and lines can be 
divided into 

1. internal lines, and vertices that connect internal lines only (blue elements in the picture); 

2. external lines, and the vertices attached to them (black elements in the picture). 

For type-1 elements, the momenta flowing into the vertices and in the propagators are all integrated over. 
Assuming we cut off the loop integrals at energies of order of, or lower than, A — which is the strong coupling 
scale of the theory — this part of the diagram contributes powers of A times the (positive) powers of A// 
suppressing the interaction vertices, if they come from expanding the exponential in ©0,2,0: or from 00,4.1- 
That is, in the parameterization of (16), it contributes 

TO = 0, n = 0, g>0. (19) 

For type-2 elements: The cubic vertex Ci,3,o always puts two derivatives on the external leg, because of the 
standard non-renormalization property [13, 15, 16], so that it contributes 

m = l, 71 = 1, q = . (20) 

The quartic vertex Ooa,i can attach to one or two external legs; it contributes 

m = 0, n = l,2, q = l. (21) 

Finally, the vertices we get by expanding the exponential in Co. 2,0 can have at most two external legs with 
one derivative each, and a generic number of external legs with no derivatives, each accompanied by a 1//. 
We can focus on the case in which there are no external legs without derivatives, since those will re-sum into 
the correct exponential structure dictated by scale invariance, which is completely fixed once we determine 
the derivative structure of the operator we generate (a particular example is the term without external 
derivatives, which will scale as ^ A'^e'^'""^^). We thus get 

m = 0, n:=l,2, q>l. (22) 

We are interested in the combination k of (17). We see that all contributions (19) — (22) have 

'?+y-2>0, (23) 

so that the whole diagram has an overall > 1, i.e., its contribution on the solution is negligible compared 
with the terms we started with. 
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Note that we are envisioning computing quantum corrections around the Poincare invariant background, 
but this background has a ghost. To justify this, we appeal to the fact that the quantum effective action does 
not depend upon which background it is computed [20]. In fact, in apphcations to spontaneous symmetry 
breaking, the true vacuum is not known a priori. Operationally, the ghost should not cause any problems in 
the computation of the effective action — the presence of a ghost in this case only changes the sign of the 
propagator and of the ie prescription (in fact, the action can be made ghost-free with a flip in overall sign, 
since we are not yet considering coupling to other matter). The infinities due to the runaway instability of 
the vacuum come from phase space factors, not from the amplitudes [21, 22], and should not be a problem 
here. 

4 Perturbations and their symmetries 

In Galilean Genesis, a scale invariant spectrum of perturbations is generated by a light spectator scalar [7], 
which we call a. Given the 5*0(4, 1) symmetry of the background solution (2), its dynamics are identical to 
what they would be in de Sitter space, where light scalars acquire a scale invariant spectrum [17, 19]. This 
spectrum can be later converted into adiabatic perturbations through a variety of mechanisms [23]. With 
our new action (1), we have broken explicitly the original 5*0(4, 2) symmetry of the model, down to Poincare 
plus dilation invariance, which does not contain the de Sitter isometry group 50(4, 1) as a subgroup. One 
might worry that our predictions for cosmological observables are impaired by the absence of full de Sitter 
symmetry. 

It is easy to see, however, that dilation invariance is all that is needed to preserve the scale invariance 
of a correlation functions'^ Indeed for a light scalar (whose lightness can be protected by an approximate 
shift symmetry), dilation invariance forces the quadratic action to take the form 



where A and are constants [24]. Notice that, in contrast to the 50 (4, 2) -invariant case, there is no reason 
to expect a to travel on the lightcone; a dilation invariant operator of the form (dirda)'^, for example, when 
evaluated on the solution (2), only affects the &^ term. The action above yields a scale invariant spectrum, 
in the same way as in models of inflation with speed of sound different from unity [25]. This is easy to see 
via symmetry arguments [7], as dilation invariance forces the 2-point function {a{t,0)a{t,x)) to depend only 
on the ratio \x\/t. As the field is massless, its wavefunction becomes time-independent at late times, and 
therefore also independent (up to logarithms) of the spatial separation. This argument extends to higher 
order correlation function of a, which will all be scale-independent as in inflation. We stress that also in 
inflation it is the dilation symmetry [rj, x) — >■ A(77, x) which is the origin of the scale invariance of correlation 
functions. Indeed, only translations, rotations and (approximately) dilations are good symmetries of the 
inflaton background^ [24] . 

While the spectrum is fixed by dilation invariance up to the overall normalization, the reduced symmetry 
of the solution will show up if we look at higher order correlation functions. In the original "luminal" Genesis, 
the action and the background solution enjoy the full 50(4, 1) group of symmetries, so that the dynamics 
of a are endowed with this symmetry, in strict analogy with what happens to a test scalar in de Sitter [24] . 
In the asymptotic future, the 50(4, 1) group acts as the conformal group in three Euclidean dimensions, 
so that the correlation functions of a are not only scale invariant but fully conformally invariant (with a 
transforming as a primary field of dimension zero, in the massless limit). In the present case, on the other 
hand, the correlation functions of a are only invariant under rotations, translations and dilations. These are 
the same symmetries of correlation functions as in single-field inflation. 

^To avoid confusion, we will call 'dilation invariance' the symmetry under rescalings 



while we will reserve 'scale invariance' for what cosmologists mean when they refer to scale-invariant correlation 
functions: that they depend at most logarithmically on distances. 

^Dilation symmetry is an approximate symmetry of inflation if all quantities, apart from the scale factor, are 
approximately constant. This "slow-roll" condition is related to an approximate shift symmetry of the inflaton. 




(25) 



a;" ^ Xx" , 



e — )■ A e 



(24) 
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These statements applied to the Unearly reahzed symmetries. But our subluminal Genesis differs both 
from the luminal case and from inflation at the level of non-linearly realized symmetries. In the original 
(luminal) model we start from an S'0(4, 2) symmetry spontaneously broken down to 50(4, 1), a special case 
of the general conformal scenario of [17]. The broken part of the symmetry is non-linearly realized on tt and 
the implications of this non-linear realization are studied in [26]. Here, on the other hand it is the broken 
part of the Poincare group only that is non-linearly realized: time translations and Lorentz boosts. This 
will give different constraints on the correlation functions involving tt. One point in common is that the two 
models both non-linearly realize time translations: for instance, as discussed in [6, 17] and as mentioned 
in Sec. 2, this fixes the mass of tt excitations. Notice that the pattern of non-linearly realized symmetries 
is also quite different from single-field infiation, where the inflaton itself non-linearly realizes the 50(4, 1) 
group [27-30] . The symmetry breaking pattern for each class of models is shown in the table below. 



Model 


Symmetry 


Linearly realized 


Single-field inflation 


50(4,l)x shift 


150(3) X dilations 


Inflation spectator 


50(4,1) 


50(4,1) 


Conformal 


50(4,2) 


50(4,1) 


Subluminal GG 


/50(3, l)x dilations 


150(3) X dilations 



Table 1: Symmetries of various cosmological models. The second line refers to a test field during infiation, 
with negligible coupling with the inflaton [24]. The third column indicates the symmetries linearly realized 
on the cosmological background: in all cases it contains rotations, translations and the dilation symmetry 
responsible for the observed scale-invariant spectrum. 

We conclude that, although all these models are degenerate at the level of the spectrum of perturbations, 
measuring higher order correlation functions would shed light on the symmetries of the system that gives 
rise to density perturbations, and would distinguish our subluminal Genesis from the luminal case and from 
inflation. 

Before closing this section, it is worth pointing out two peculiarities of our symmetry breaking pattern. 
First, the de Sitter solution (2) is more symmetric than it should be, in the sense that it has some accidental 
symmetries that are not present in the action (1). Indeed, we have been calling it 'de Sitter solution' precisely 
because it is invariant under the de Sitter isometry group, 50(4,1). This, however, is not a subgroup of 
the symmetry group that leaves the action (1) invariant, which is just the Poincare group supplemented 
with dilations. The dynamics of tt perturbations and of other flclds (cr) that couple to tt are invariant under 
the symmetries of the action, with the spontaneously broken ones realized non-linearly. Possible accidental 
symmetries of the solution we expand about do not translate into accidental symmetries for the dynamics 
of TT and a excitations — hence, the analysis in this Section is unaffected by this subtlety^. 

Second, the fact that the de Sitter solution (2) is invariant under dilations, 

7i'ds(a;) 7rds(Aa;) + log A = 7rds(2:) , (26) 

means that dilations are realized linearly on tt's perturbations, 

(f{x) = Tr{x)~TTds{x)^f(\x), (27) 

even though they are realized non-linearly on the original tt held. This is quite uncommon. Usually, when we 
expand an action about non-trivial held configurations, we break some of the original symmetries, and for 
those broken symmetries we go from linear representations to non-linear ones. Here, instead, we start with 
Poincare transformations realized linearly and dilations realized non-linearly on tt, and we end up with the 
opposite. The Poincare invariant solution = 1 breaks dilations; the non-trivial solution we consider, (2), 
breaks some of the Poincare generators (Lorentz boosts, time translations), but restores dilations. 

'^Notice that having enhanced 'background' symmetries is not unusual. A simpler example is given by an 
anisotropic model with Lagrangian L = (p — ci{dx<t>)^ — Cy{dy4>)'^ — c1{dz4>)^ ■ The background solution 4> = const, is 
invariant under the full Poincare group, even though the action is not. 
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5 Junction with standard cosmology 



The Genesis phase must eventuaUy be foUowed by the standard radiation dominated phase. We wih not 
provide an exphcit model for the mechanism that reheats the universe (see [31]). Rather, we will focus on 
general energetic considerations and highlight a peculiarity of our scenario both in its luminal and subluminal 
version. 

Regardless of the details of how the galileon energy gets converted into radiation, one would naively 
expect that for a sufficiently fast transition — much faster than the Hubble rate — the energy density and 
therefore H are continuous across the transition. The logic is that the conservation of the stress-energy 
tensor gives 

p=~'SHip + p), (28) 

so that the overall variation of p is small if the transition time is short compared to H^^. Note that this 
argument implicitly assumes that p is not parametrically larger than p, which is usually the case. In our 
Genesis phase, however, we have \p\ 3> p, so that the argument does not go through. This is the same as 
saying that the rate of variation of p is not set by H, as in a standard phase, but by 1/t, which is much 
larger. 

One might be led to conclude that only a transition that is rapid compared to t guarantees the continuity 
of H. It turns out that this conclusion is also too hasty. The point is that, as we make the transition faster 
and faster, the pressure blows up. Indeed, from the explicit form of p and p given in (6) and (7), we see that 
the last term in the expression for p, 

blows up as we make an instantaneous transition in which tt changes by a finite amount. To be concrete, let 
us imagine that when tt reaches a certain value tt* , there is a sudden upward change in the potential energy 

T/(7r) = 6'(7r-7r,) . (30) 

In this way, (part of) the Galileon energy is transformed into potential energy, which later can be easily 
converted into radiation or anything we like. In the absence of the singular contribution to the pressure, 
p would be continuous across the transition, and in principle one could dial the value of Vq to get all the 
galileon energy converted into potential energy. However, it is easy to check that the equation of motion for 
TT is (obviously) equivalent to (28) and gives 

P = -3H ■ {psing + finite) , (31) 

where 'finite' stands for all contributions to p + p other than Psing, which are finite — though potentially 
ike V{tt) — at the transition time. Using the I 

be rewritten as 



discontinuous, like V{Tr) — at the transition time. Using the Friedmann equation, = p/3AIpi, this can 



37(^-7T7^^i2T^) ^fiiiite. (32) 



dt\ 9{l + a) H'^ M^i^ 

If we now integrate this expression in time from slightly before to slightly after the transition, we see 
that it is the quantity 

2 P TT^ 

9{l + a)HlMl ^^^^ 

that remains constant — rather than p or H — across an instantaneous transition. In Appendix B, we 
show how this result can be obtained covariantly through Israel-like junction conditions. Note that p = 
on our Galilean Genesis solution (2), hence _ff = to zeroth order in 1/A/pi. To obtain the leading non-zero 
contribution to H{t), we can integrate HAIp^ = —{p + p)/2 to obtain 

In the NEC-violating range (9), this describes an expanding universe from an asymptotically static state, 
as it should. Substituting (2) and the above expression for H{t), we find that (33) equals ^Epf -ff- Hence, 
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assuming tt vanishes after the transition, the Hubble parameters before and after the transition arc related 

by 

I + a 

Haiter = -^before (assuming IT — > instantaneously) . (35) 

3 + a 

In principle there is no reason to restrict to an instantaneous transition. However, at any given time t, 
there is a time \t\ before H blows up, so the transition better conclude within a time at most of order t. For 
a smooth transition with this typical time scale, all terms in p are of comparable sizes, 

/2 1 

P-^^, (36) 

and p undergoes an overall variation of order one, 

Ap'^tHp^ p . (37) 

This implies 

-ffaftcr ^ ^^bcforo (assuming n ^ in At ^ t) . (38) 

Apart from a model-dependent numerical coefficient of order one relating -ffboforc a-nd i?aftcrj these results 
do not change the qualitative conclusion that our NEC-violating phase can take an initially vanishing or 
negligible Hubble rate, increase it by many orders of magnitude, and then make a transition to a standard, 
NEC-obeying radiation-dominated cosmology. However, they do offer an unexpected twist for the whole 
start-the-universe-via-NEC-violation program: not only does one have to come up with consistent NEC- 
violating mechanisms that create energy out of nothing, but one must also make sure that the bulk of this 
energy can later be passed to more standard forms of matter or radiation. 



6 The other Gahleons and couphng to gravity 

The starting point of the previous Sections was the minimal conformal Galileon Lagrangian, which includes 
only the kinetic term and the cubic interaction together with their conformal completion. The existence of 
a stable NEC-violating solution in this case forces the kinetic term to have the wrong sign [6, 7], and this 
in turn implies that the fluctuations around the Poincare invariant solution n — are ghost-like. However, 
there are five possible conformal Galileon terms [5], and one may wonder whether the presence of higher 
order interactions or their dilation invariant deformations can give rise to a stable Poincare invariant vacuum, 
without spoiling the good properties of the de Sitter solution discussed above. 

The inclusion of the higher order Galileon terms introduces an ambiguity in the choice of the action — 
and as a consequence in the definition of the NEC itself — when the coupling of the scalar to gravity is 
considered. The fact that non-minimal couplings to the metric can give rise to different stress-energy tensors 
is not surprising, the most familiar example being the "improved" energy-momentum tensor obtained by 
adding the conformal coupling —j^Rip^ to the minimally coupled, massless Xip'^ theory. What is new in the 
presence of higher order Galileon interactions is that even the definition of minimal coupling is ambiguous [32] . 
The reason is easy to understand: consider the quartic Galileon in flat space, schematically of the form 
(971)2(9^71)2, then commute some derivatives and integrate by parts to obtain the same structure with a 
different Lorentz contraction. Now let us minimally couple the two (equivalent) structures by promoting 
partial derivatives d to covariant ones V and contracting all indices with g^^^. We can go back from the second 
structure to the first one by integrating by parts again but now, since the theory has higher derivatives, there 
can be a non-trivial commutator [V, V] which is proportional to the Riemann tensor, thus a term of the form 
i?(97r)^ can appear*: the two "minimally coupled" interactions are not equivalent and they give different 
contributions to T^u- We will discuss below which one among all possible minimal couplings is the most 
convenient for our purposes. 

So far we have concentrated on the case(s) of minimal coupling, however it has been argued that this 
choice may not be the healthiest. Indeed, in the presence of dynamical gravity, the contribution to the 

*In the case of the cubic term this ambiguity is absent because this term does not have enough derivatives to 
generate a non-trivial commutator. 
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equations of motion for the scalar and for the metric perturbations given by the quartic and quintic GaUleon 
do not remain of second order, but also contain three derivatives for any of the possible choices of minimal 
coupling [lO]''. Operators with three derivatives are not obviously associated with extra ghost-like states — 
at worst, one needs to impose one additional initial condition, which sounds like half an additional degree 
of freedom. Still, it is hard to believe that they do not impair the validity of the theory in situations where 
they become relevant. One can avoid them altogether by adding suitable non- minimal couplings [10]. They 
are truly non-minimal in the sense that they do not correspond to any of the possible minimal choices 
discussed above, indeed one needs terms that cannot be generated by commuting covariant derivatives. 
These non-minimal couplings can give a different stress-energy tensor. However, the calculation of T^i/ is 
considerably more complicated, even restricting to the homogeneous and isotropic case, as one cannot use 
the tricks discussed in [6]. Even these non-minimal terms are not unique, and there are different choices 
which yield second order equations and reduce to the Galileons on flat space [33]. A preferred choice for the 
non-minimal terms can be derived using the brane construction of [34-36]. We further explore this choice 
and the associated stress tensors in Appendix A. 

From an effective field theory (EFT) standpoint, the three-derivative terms are not pathological insofar 
as they can be treated as a small perturbation. The cutoff of the EFT must be lower than the scale where 
these higher-derivative terms become important. This scale is easy to estimate for our cosmological solution. 
The Lagrangian with one of the minimal couplings to gravity takes the form [6] 



where (j) = Hqc^ , so that the quadratic action for the perturbations, 



^'-ffS.^V|Mj,«, (3!,) 



7 + -0 , 9tJ.u = VtJ.u + hf,y , (40) 



t 

reads schematically 

^ [{d^l^r + d^i^dh] + MUdhf , (41) 

where we neglected mixing terms with fewer derivatives, to be discussed shortly. Once we choose canonical 
normalization for the diagonal terms, the mixing takes the form 

TFT^^'V'cSftc . (42) 

It becomes as important as the diagonal terms at an energy scale 

-Einax = —p-Ho , (43) 

which sets an upper limit on the UV cutoff of the EFT^*^. Notice that this scale depends on Mpi, so that it 
is parametrically different — though not obviously smaller or larger — than the other UV scales discussed 
in [6], such as the strong coupling scale A ^ (/iJp)^/^. The background solution has a typical energy scale 



^Again, since [7] and the present paper focus on the simplest scenario with only the cubic conformal Galileon, 
their analyses are not affected by the following discussion. 

^°The effect of the three-derivative terms is not always small as it happens on this background. We can repeat the 
analysis of the previous paragraphs in the case where the Galileon field is responsible for the present acceleration [5] 
and concentrate on the spherical solution sourced by a localized object of mass such as a star. Focusing on the 
region outside the Vainshtein radius, the quadratic action for the perturbations 

, . M, 1 

■n = TH)(r) + if = ifi , ff^f = '7m-' + "mi- , (44) 



is schematically (taking / ~ Mpi) 



MUdipf + -^{duofdhd'^ip + Mliidhf 



(45) 
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d(j)/(j) ^ 1/t which, neglectmg gravity momentarily, sets the freeze-out frequency of scalar perturbations [7]. 
Therefore cosmological observables can be consistently calculated within the EFT if this scale is smaller than 
-Emaxj which is the case for sufficiently early (negative) times: 

1^1 » to ^ -^H^^ ■ (47) 

This is the same regime in which we can neglect the effect of gravity on the background solution [7] . Indeed 
going to early times is equivalent to sending Mpi to infinity, which makes all the gravitational effects — 
including the mixing — weaker and weaker. In particular, in this regime the freeze-out of perturbations is 
dominated by the (j) background, as assumed above. It is easy to check that also the two-derivative mixings 
are negligible in the same limit, scaling as t^/t ■ dtpcdhc- In conclusion the minimal coupling is perfectly 
consistent as long as we stick to sufficiently early times, in the sense of (47). If the transition to the standard 
cosmology occurs before the model is consistent throughout. In fact, we could even envision that reaching 
t ^ to and consequently probing energies of order -Bmax is what triggers this transition: new physics must be 
present at -Emax to save the consistency of the theory, and this new physics may be responsible for draining 
energy out of the Galileon sector and reheating the Universe. 

The conclusion of this digression is that as far as we are interested in our early-universe scenario, the 
minimal couplings do not give rise to instabilities below the UV cutoff of the EFT. Notice that all these 
choices (minimal couplings, "covariant improvements" [10], etc.) are also technically natural. As shown 
in [37], starting with minimally coupled Galileon terms, we generate operators of the form (9^7r)", and 
terms where two derivatives are replaced by a Riemann tensor, which are subleading. We can then use the 
minimally coupled Lagrangian that is obtained from a specific form of the conformal Galileon operators, the 
one built starting from the curvature invariants involving the conformally flat fictitious metric e^'^rj^i, [5]. 
The computation of T^^ in this case drastically simplifies [6] . 

We can finally answer the question raised at the beginning of this Section. When all the Galileon operators 
are present, the sign of the kinetic term can be reversed to give a stable tt = solution while preserving the 
existence and stability of the NEC-violating one. However, it inevitably exhibits superluminality. Around 
the Poincarc invariant solution fluctuations are luminal, independently of the breaking of special conformal 
transformation; the leading correction to the speed of propagation around weak-field deformations is given 
by (97r)^n7r, and it is superluminal [6]. The only possibility to avoid superluminality would be to set the 
coefficient of the cubic interaction to zero, but this is not compatible with the existence of a stable NEC- 
violating solution as manifested by the conditions found in [6]. This conclusion holds for a = or small 
enough. One could hope that for a large value of a, not only one has subluminality around the NEC violating 
solution, but also the Minkowski vacuum becomes healthy. Unfortunately this is not the case, as we prove 
in Appendix C. 



7 Conclusions 

Subluminal Galilean Genesis is a consistent model of early cosmology that is alternative to standard inflation. 
The time-dependent NEC-violating solution obeys all the basic theoretical consistency requirements: it is 
stable, excitations around it are largely subluminal and it comes from an action whose form is protected 
by approximate symmetries. The same symmetries lead to the production of a scale invariant spectrum of 
perturbations. The predictions of the model are not completely degenerate with the ones of inflation and a 
distinction between the two is possible at the level of higher order correlation functions. 

Various questions remain open. In our framework we could not build a theory that, besides the cosmo- 
logical solution we are interested in, describes also an healthy Minkowski vacuum. It is not clear whether this 



Going to canonical normalization one finds 



At a distance from the source of the order of the Vainshtein radius rv ~ (^^)^'''^A~^ the cutoff is comparable to the 
inverse Schwarzschild radius, -Emax ~ l/^s- This may not be problematic for the solar system, but it further limits 
the validity of the EFT. 
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is accidental or points to a subtle inconsistency implied by the violation of the NEC. To explore this issue 
one could consider different couplings with gravity and/or extensions of the Galilean symmetry. Another 
question is whether there are other theoretical consistency checks — analogous to the usual analyticity prop- 
erties of the S-matrix in Minkowski — that can be applied to a time dependent solution. Finally it remains 
open the general question of whether some fundamental properties of UV-complete relativistic theories forbid 
NEC violations, even if consistent EFTs can be written down. 
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A Non-minimal couplings and the higher Galileons 

Here we calculate explicitly the possible higher Galileons along with the set of non-minimal couplings pre- 
ferred by brane constructions. We then calculate the energy density and pressure around a dS like solution, 
and obtain a curious result showing that the NEC is violated if and only if the coefficient of the quadratic 
term around flat space is ghost-like. 

In [38] , it is shown how to add non-minimal terms to any kind of Galileon in order to preserve the second 
order equations of motion. However, the procedure is not unique - there are many different non-minimal 
terms one can add which preserve second order equations and reduce to the desired flat space Galileon as 
the metric goes to Minkowski. Different choices can give different energy densities and/or pressures, even 
around flat space. 

A way to naturally single out a particular choice of the possible non-minimal terms of [38] is via the brane 
construction of [34-36]. When constructed in this way, the Galileons have a 5-dimensional interpretation: 
the scalar, in a certain limit, is the brane bending mode of a flat brane (with higher order world-volume 
curvature corrections) probing AdSs with the zero mode of the bulk metric turned on. With the metric turned 
off, the non-linear conformal symmetry descends from the isometry group of AdSs . Following through this 
construction, we find the following set of non-minimal couplings for the five conformal Galileons, 



C5 



47r 



-4tt/ 



ge-^-(97r)^ -[H]^ + [H 



{dnf [Uf - 3[H][h2] + 2[h3] - 3(97r)2([n]2 - [H^]) 



30 3 



1. 



1 



(48) 



Some explanation of the notation is in order; tli6 metric is the associated covariant derivative is V^, and 
G^u = R^i, — \Rgfiv is the Einstein tensor. Meanwhile, H is the matrix of second derivatives H^^ = V^tt. 
For traces of the powers of H we write [H"] ee Tr(H"), e.g. [H] = V^V^tt, [H^] = V^Vi.TrV^VTr, where 
all indices are raised with g^". We define the contractions of the powers of H with Vtt using the notation 
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[tt"] EE Vtt • n"~2 . Vtt, e.g., [tt^] = V^TrVTr, [tt^] = V^ttV V^ttV^tt, where all indices are raised with 
respect to g^^'^ . 

The only two terms used in the present paper and in [7] are £2 and £3, and these have no non-minimal 
terms or ambiguities in the minimal coupling. 

Varying the Lagrangians (48) with respect to the metric and then specializing to time-dependent con- 
figurations TT = iT{t) on flat space g^j/ = rj^^, we find the following energy densities and pressures, 



P2^\e'-i^\ P2 = \e'^i^\ (50) 

3 1 

Pz = ^Tr" , p3 = -TT^ (tt^ - 4^) , (51) 

P, = \e-^-i^\ P4 = ^e-2-^4(^2_2^) ^ (52) 

P5 = le-'^'ir^ , P, = (5^2 _ . (53) 

For the dS configurations of interest 7r(t) = — log(— iJoO; given the general lagrangian C = Y^i=i'^i^i 
with arbitrary coefficients Ci , the equation of motion for flat space reduces to 

ci + 2c2Hl + ScsH^ + 3ciHl + "^^^ = Q , (54) 

so a solution exists for any positive real valued root to this polynomial in Hq. 
The total energy density and pressure on this solution are 

El 13 S-ffn SiJn , , 

El 13 ^i/n ^i/n , , 

^C.P. = -Cl:|^+C2^-C3^-C4^-C5^. (56) 



1=1 



Remarkably, all the terms except those from the quadratic part £2 cancel in the combination p + p, 

P + P=^, (57) 

and this is negative if any only if C2 is, 

NEC violation 4^p + j3<0^C2<0. (58) 

Thus, no matter what higher order Galilcons are included, violating the NEC on a genesis-like solution 
requires a wrong sign kinetic term around the Minkowski vacuum. Note that this cancellation for the higher 
order terms is something special about the special non-minimal couplings in (48) coming from the higher- 
dimensional brane construction, and would not be true for minimal couplings or other possible non-minimal 
couplings. In particular, it is not true of the choices made in [6]. 



B Covariant Junction Conditions 



In this Appendix we derive from a covariant standpoint the junction condition to standard cosmology ob- 
tained in (33). In the limit of an instantaneous transition, we can approximate the transition event as a 
space-like surface located at some t — t^,. The action for i < is given by 

-g + fe'^dnf + ^.{dnfUn + ^(1 + a){dnf'\ 



S[t < t,] = 



t<t. 



2A3" 



d^xVh 



(59) 
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where hij denotes the induced metric on the t = t^, boundary, Kij = Cnhij /2 is its extrinsic curvature, and 
Cn denotes as usual the Lie derivative with respect to the unit time-Hke vector normal to the boundary. 
The boundary action includes the Gibbons-Hawking term and its galileon cousin [39-41] necessary to have 
a well-defined variational principle. Similarly for t > t^,. 

Stationarity with respect to variations of the metric gives Einstein's equations, together with an Israel 
junction condition at < = i*: 



A 



= 0. 



(60) 



In other words, the quantity in square brackets is conserved across the transition. Specializing to cosmology, 
if* = HS^j and tt = 7r(t), and using (3), this reduces to 



A H- 



2 P 
9(1 + a) Hi 



= 0, 



(61) 



which reproduces (33). 



C Healthy Minkowski vacuum and higher Gahleons 

In this Appendix we show that even considering all the four conformal Galilean operators (with minimal 
coupling with gravity) deformed by the addition of a term (Stt)^ proportional to a, as in eq. (1), it is 
not possible to have a NEC violating solution and, at the same time, a completely healthy theory around 
Minkowski. 

We make the following requirements on the theory: 

1. ci = 0, to allow for a Minkowski solution. 

2. C2 > to avoid ghosts around Minkowski. 

3. C3 — 0, since the DGP term, (97r)^n7r, always induces supcrluminality around weak field solutions in 
Minkowski, independently of the sign of the operator [6]. 

4. a > 0, both to close the lightcone of perturbations around the NEC violating solution and to have an 
healthy 2-to-2 S-matrix in Minkowski [6]. 

5. C4 < to avoid supcrluminality around weak field solutions in Minkowski. Indeed, after setting C3 — 0, 
the leading corrections to the lightcone come from the conformal breaking (Stt)"* and from the quartic 
Galileon. Comparing the coefficients of the two operators (with a ^ 1), it is easy to see that the 
quartic Galileon dominates when the weak field solution is characterized by a length scale which is 
much smaller than Hq, d Hq. Focussing on this regime, we consider the linear equations of motion 
around a classical background ttq and we are interested in the modification of the lightcone for the 
perturbations. The quartic Galileon gives the following contribution to the equation of motion: 

{anf - 3n^(9^a,7r)2 + 2id,,d,T:f . (62) 

The first term does not change the lightcone aperture. As the linear equation of motion is Ottq = 0, 
in the second term we are forced to put the two d^duir legs on the background, so that also this term 
does not change the speed of the fluctuations. The third term evaluated on a static background gives 

(didjTTndjdkTTo) Okdiir , (63) 

and the matrix in brackets is positive definite. This implies that the coefficient of this operator must 
be negative to give subluminal propagation of perturbations. 

Let us see whether these conditions are compatible with the existence of a NEC violating solution. The 
equation of motion for a 'de Sitter' solution can be easily derived from [5] using the following trick. For 
a time dependent solution the term (97r)^n7r does not contribute to the equations of motion, so that the 
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results of [5] can be used simply replacing the coefficient of the cubic conformal Galileon C3 with the one of 
{dirf ("). This gives 

- 2c2 +3aff^ -3c4iJ^ + ^csi?^ = . (64) 

Also the arguments of [6] to calculate the stress-energy tensor can be used; indeed the cubic Galileon is 
anyway treated separately. We get 

where in the last step we used the equations of motion above. The quantity in brackets should be positive to 
violate the NEC, but all the terms are constrained to be negative because of the conditions discussed above. 
This shows the impossibility to get an healthy Minkowski background, even before imposing the additional 
constraint from the stability of the NEC violating solution we found. 

The previous discussion brings up a possible loophole in the arguments of the paper. To have an 
healthy NEC-violating solution, we considered the breaking of 5*0(4, 2). But would it be possible to achieve 
the same result, preserving the original symmetry, by simply requiring that the coefficient of the cubic 
Galileon vanishes around the NEC-violating solution and that the quartic Galileon has the healthy negative 
sign discussed above? Following the recipe of [5, 6] it is straightforward^^ to check that this is indeed 
possible (though not compatible, as in the rest of the paper, with an healthy Minkowski vacuum). However, 
one quickly realizes that the situation is not completely satisfactory. A negative coefficient of the quartic 
Galileon ensures the subluminality of perturbations around time-independent solutions. If this may be a 
good criterion around Minkowski with non-relativistic sources, it is clearly not enough for the NEC- violating 
time-dependent solution, where space- and time-dependent perturbations must be considered. The last term 
of (62), however, induces also superluminal corrections to the lightcone once time-dependent backgrounds 
are considered. 
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